Cooper pair insulator in amorphous films induced by nanometer-scale thickness variations 
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Unusual transport properties of superconducting (SC) 
materials, such as the under doped cuprates|l|, low di- 
mensional superconductors in strong magnetic fields 1 2 1, 
and insulating films near the Insulator Superconductor 
Transition (IST)|3|, have been attributed to the forma- 
tion of inhomogeneous phases. Difficulty correlating the 
behaviors with observations of the inhomogeneities make 
these connections uncertain. Of primary interest here 
are proposals that insulating films near the 1ST, which 
show an activated resistance and giant positive magne- 
toresistance, contain islands of Cooper Pairs (CPs)(!H2l- 
Here we present evidence that these types of inhomo- 
geneities are essential to such an insulating phase in amor- 
phous Bi (a-Bi) films deposited on substrates patterned 
with nanometer-sized holes. The patterning induces film 
thickness variations, and corresponding coupling constant 
variations, that transform the composition of the insula- 
tor from localized electrons to CPs. Analyses near the 
thickness-tuned ISTs of films on nine different substrates 
show that weak links between SC islands dominate the 
transport. In particular, the ISTs all occur when the link 
resistance approaches the resistance quantum for pairs. 
These observations lead to a detailed picture of CPs local- 
ized by spatial variations of the superconducting coupling 
constant. 

A central issue in the field of insulator to superconduc- 
tor transitions concerns the nature of the insulating phase. 
In most homogeneously disordered thin films, the insulator 
is composed of unpaired electrons that are weakly localized 
by disorder (SHIIl. A less familiar state which supports rem- 
nants of superconductivity appears in other homogeneously 
disordered films, such as In Oxide and TiN lfT2l [T3ll . This 
so-called Cooper Pair Insulator (CPI) state exhibits activated 
transport and a giant positive magnetoresistance, which sug- 
gest that localized Cooper pairs participate in transport. How 
this bosonic insulating state develops has not been deter- 
mined. Theoretical descriptions of the CPI phase often pre- 
sume that disorder or magnetic field creates islands of Cooper 
Pairs (CPs) embedded in a non-SC backgroundfTl [T4lfT6ll . In- 
deed, scanning tunneling spectroscopy experiments show evi- 
dence of incipient puddling in some films close to the IST lfTTl . 
Since Cooper pairs are not evident in most other insulating ho- 
mogeneously disordered thin films, the factors responsible for 
the change in behavior are generally unclear. 

Recently, it was established that amorphous Bi films can 
be induced to form a CPI state instead of the unpaired insu- 
lator state by depositing them on an anodized Al oxide sub- 
strate patterned with a nanohoneycomb (NHC) array of holes 



(see Fig. [TJa) and (b) insets). The holey substrates have non- 
uniform topography that causes thickness and, consequently, 
superconducting coupling constant variations throughout the 
film. The insulating NHC a-Bi films near the 1ST exhibit the 
characteristic transport signatures of the CPI phase (18] QjO 
. Moreover, the insulators are confirmed to contain Cooper 
pairs through Little-Parks like magnetoresistance oscillations 
induced by the hole array |[T8l . These films provide a unique 
example of a system developing a CPI phase as a result of the 
introduction of inhomogeneities. Establishing that well de- 
fined islands form at the 1ST will put important constraints on 
models of the Cooper Pair Insulator phase. 

Here we focus on the systematics of thickness tuned ISTs 
in these films. We pay particular attention to how the critical 
thicknesses for the onset of conduction (doc) an d supercon- 
duction (c?ist) var Y with substrate qualities and also point out 
the relative invariance of the critical normal state sheet resis- 
tance C^ist)- The ISTs, shown in Fig. [TJa) and (b), differ 
significantly from those of their unpatterned, or "reference," 
counterparts (see Ref. for an example). The R(T) evolve 
with deposited thickness, d dep , from a monotonic rise as T — » 
that is exponential rather than logarithmic, to develop a reen- 
trant dip not seen in reference films, to finally drop steadily 
into the SC state. The exponential rise indicates thermally ac- 
tivated transport, R(T) ex e T °/ T (\M- For all the ISTs, the 
normal state sheet conductance, G(SK) = 1/R(SK) grows 
linearly with deposited film thickness, d dep , starting at a crit- 
ical thickness, <iQc(Fig. [TJc)). This linear behavior reflects 
homogeneous rather than granular film growth, which would 
appear as an exponential dependence! 20]. Consistently, the 
slope, dG(SK)/dd dep , is smaller and the intercept, do c , is 
larger for the NHC films compared to reference films. Tq de- 
creases continuously on approaching the 1ST (also shown in 
Refs. (mdUEEl). We define the critical thickness for the 
1ST, disT> as the thickness at which T linearly extrapolates 
to zero. Notice that similar to d& c , d\^T is larger for NHC 
films than reference films. Finally, Rist shows little varia- 
tion over nine different substrates that exhibit a wide range of 
<4st (see Fig. fTJe)). The i^isTS cluster around Rist ~ 3Rq, 
where Rq = h/(2e) 2 , the quantum of resistance for pairs, is 
the critical sheet resistance for SC in reference films. 

Topographical analyses of AAO substrates using AFM im- 
ages show that the minimum film thickness in critical NHC 
films (at d dep = doc and c?ist) correlates with critical thick- 
nesses in reference films, d T Q C and d[g T . As pointed out 
previously [22], the NHC substrate topography h(x, y) shows 
regular variations around holes, which makes the thermally 
evaporated film thickness a function of position, d(x, y). For 



2 



G 

CD 
O 

c 



CO 


DC 

"55 

_c 

CO 





3 °n .Q-- 



d dep (nm) 



- 1 



— i— r- 1— i— i- 

. e 


i ■ 










. 1 1 . . 


1 ■ 


1 ■ ■ ■ 



J3 



2 JJ 

1 ° 



0.1 1 10 0.5 

Temperature (K) 



1 1.5 

d IST (nm) 



Figure 1 : Characteristics of thin films on nanohoneycomb sub- 
strates (a) and (b) Insulator Superconductor Transitions tuned by 
thickness of a-Bi on more and less ordered hole arrays. Insets: post- 
experiment SEM images. Center-to-center hole spacing: 100 ± 5 
nm. Radii: 27 and 28 ± 3 nm. The scale bar spans 200 nm. (c) Nor- 
malized conductance (Coo = l/81kQ) and (d) activation energies 
for films on four different NHC substrates. Similar symbols repre- 
sent data from experiments on side-by- side NHC substrates. Unfilled 
and filled symbols represent data from more and less ordered hole ar- 
rays, respectively. Data from the experiment in (a) and (b) are shown 
as circles. For squares, rhoie = 16 (unfilled) and 15 (filled) ± 3nm. 
Reference film data are shown as grey triangles in (c). The dotted 
lines are least- squares linear fits to the data, and the arrows indicate 
critical thicknesses for conductivity and SC. Open arrows indicate 
critical thicknesses on a reference film: 0.5 nm for conductivity and 
0.7 nm for SC. (e) Normal state (T = 8K) critical sheet resistance 
for the last insulating and first SC films of the 1ST of nine different 
NHC substrates. The horizontal dashed line at 3Rq is a guide to the 
eye. 



a uniform flux of atoms impinging parallel to the average nor- 
mal to the substrate these are given by: 



d(x, y) = d dep 



(1) 



Fig. |2|a) shows the topography of the A AO substrate. Six 
peaks, visible as bright spots, surround each of the holes. 
A line scan between two neighboring peaks shows a single 
central minimum in the height (Fig. |2jb)). Using Eqn. [T] 
we calculate thickness profiles along this link corresponding 
to the onset of conduction (d dep = c?Gc)» an insulating film 
(d Gc < d dep < d IST , bold in Fig. [lja)), and the onset of SC 
(d dep = disT) (see Figs. [5Jc-e)). All show two minima that 
correspond to the steepest regions of the height profile. At the 
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Figure 2: NHC film topography (a) AFM image showing the topog- 
raphy of the substrate in Fig. [TJ a )- The scale bar (white line) spans 
100 nm. (b) Height profile along the trace (black line) in (a), (c)-(e) 
Film thicknesses along the trace in (a) for d dep at (c) the onset of 
conduction, (d) a representative insulating film (bold in Fig. [TJa)), 
and (e) the 1ST. Pink and blue lines indicate the critical thicknesses 
for conductivity and SC, respectively, for a reference film, (f) Aver- 
age link thickness for NHC films at their transitions to conductivity 
(dGc, pink) and SC (cZist, blue) for six different substrates. Critical 
thicknesses for SC are shown by plotting the link thicknesses of the 
last insulating and first SC films spanning the 1ST. 



critical thicknesses, do c and g?ist» the minima coincide with 
critical thicknesses in reference films, d T Q c (pink line at 0.5 
nm) and d\f T (blue line at 0.7 nm). These observations inti- 
mate that NHC films conduct or SC when the thinnest regions, 
the weak links, reach the appropriate critical thickness. Cal- 
culations of the weak link thicknesses averaged over each of 
the six substrates analyzed, (? mks , strongly support this sug- 
gestion. Fig. [2^f) compares d 1111 ^ at NHC film critical thick- 
nesses with the critical thicknesses for conductivity and SC 
of reference films. It is evident that for all but one of the sub- 
strates d} inks at critical d dep corresponds to critical thicknesses 
in reference films. 

Full two-dimensional maps of the local film thickness, 
d(x, y) reveal a compelling picture: SC and normal regions 
coexist in NHC films through the 1ST (Fig. 3). These maps 
were generated by applying Eqn. [T]to the AFM image in Fig. 
|2ja) for the series of depositions in Fig. [2jc), (d), and (e). The 
color scale indicates regions that are SC, resistive, and insu- 
lating according to their behavior at those thicknesses in a ref- 
erence film. It is evident that SC islands form and grow where 
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Figure 3: Film thickness variations in NHC films Calculated thickness variations in films deposited on the NHC substrate of Fig. [TJa) and 
Fig. [2]a). Film thickness colors represent behavior on a reference film: white is non-conductive (d < cIqI), pink is conductive (d > d T Q f c ), and 
blue superconduct (d > d\s T ). (a)-(c) Film thickness variation plots (a) at the onset of conduction, (b) for a representative insulating film, bold 
in Fig. [l} a X and (c) at the 1ST. The critical temperatures for SC in reference films (T c ) of these thicknesses are shown on the right. T c = 
marks the 1ST and the maximum T c seen in each of the panels is (a) 0.4K, (b) 1.3K, and (c) 2.2K. (d) T c (in Kelvin) contours shown for a 
section of the film at d dep = disT- (e) Illustration of a wire array of resistors, R hnk , connecting SC islands. The black segments show the 
basic element of the array. 



the substrate is flattest, such that there are generally 12 islands 
surrounding each hole, one for each peak and valley. Through 
the 1ST, the SC islands are larger than the coherence length 
in these films, £ ~ 15 nm. The right axis of the figure further 
converts d(x, y) to a local T c map, T c (x, y), using the relation- 
ship between T c ref and d ref observed here and in many amor- 
phous filmsystemsilEl: T c ref = T c bulk (l - dg T /d dep ). For 
a-Bi, T h c ulk = 6 K and d\% = 0.7 nm[24j. Note that the 
T c maps are better thought of as coupling constant maps since 
they do not include proximity effects. Each of the islands has 
a larger T c at its center than on its edges and resistive film 
regions surround the islands. With increasing thickness, the 
intervening resistive regions shrink until the SC islands coa- 
lesce (see Supplemental Information). The T c contour plot, 
Fig. |3jd), emphasizes that T c in the regions between islands 
at their coalescence is finite, but much less than the maximum 
at the islands' centers. Thus, these regions likely act as weak 
links. 

The color maps suggest a network model for NHC films 
that provides a clue to the origin of the apparently universal 
normal state resistance at the 1ST. In the model, the SC is- 
lands become highly conductive regions connected by resis- 
tive weak links as shown in Fig. [3je). The resistance of the 
basic building block of this network, which is equivalent to 
the macroscopic sheet resistance, is 3R hnk . Thus, according 
to Fig. [IJf), at the 1ST « R Q . Interestingly, theory E) 
and experiment! 26 ] maintain that arrays of small, resistively 



shunted Josephson Junctions in the extreme quantum limit un- 
dergo an 1ST for R^ w Rq . Applying these dissipation driven 
1ST models to NHC films seems plausible as the SC islands 
are small enough to be in the extreme quantum limit and the 
resistive regions surrounding the islands can act as the dissi- 
pative elements. 

We have shown that the CPI phase in NHC a-Bi films is 
composed of superconducting islands. These islands result 
from variations in the film thickness on scales larger than £ 
that lead to inhomogeneities in the superconducting coupling 
constant. Our observations support conjectures that CP is- 
landing gives rise to the activated resistance and giant mag- 
netoresistance observed in this and other thin film systems, 
such as InOx and TiN. For these structure-free films, how- 
ever, an alternative mechanism for island formation is re- 
quired. Whether this mechanism is spontaneous electron- 
electron interaction driven phase separation l27l . the fractal 
nature of the electronic wavefunctions close to the metal to 
insulator transition! 6], a random Gaussian potential! 14], un- 
derlying structural inhomogeneities in the material I51I281L or 
yet to be determined is an open question. We anticipate that 
the detailed picture of the superconducting inhomogeneities 
presented here will aid microscopic calculations of these and 
other characteristics, such as the tunneling density of states, 
of the CPI state. 

Methods: 

NHC a-Bi thin films are grown by thermally evaporating Sb 
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and then Bi onto an Anodized Aluminum Oxide (AAO) [29] 
substrate held at 8K in a dilution refrigerator cryostat and 
measured in situ. Reference films are grown simultaneously 
on a fire-polished glass substrate. Film thicknesses are mea- 
sured with a quartz crystal microbalance. The semiconduct- 
ing Sb underlay er (d ~1 nm) promotes the surface wetting 
of the subsequently deposited Bi. Previous STM analyses 
of quench condensed metal films on semiconductors indicate 
that the metal atoms stick where they land to form a homoge- 
neous amorphous film l30l . Thus, the films assume the holey 
geometry of the underlying AAO substrate (see inset Fig. [TJa) 
and (b)). In the nine different AAO substrates used in these 
experiments the average hole radii range from 14 to 28 ± 3 
nm with average center-to-center hole spacing of 100 ± 5 nm 
in each substrate. The order in the hole arrangement was in- 
tentionally varied using two different anodization processes 
(see insets Fig. [TJa) vs. (b)). Sheet resistances as a function 
of temperature, R(T), were measured on 1.5 mm square ar- 
eas defined by predeposited Ge/Au contacts using four point 
AC and DC methods. The geometry of each of the nine hole 
arrays was characterized using Scanning Electron Microscope 
images and the surface topography of six of the nine was mea- 
sured using tapping mode AFM over a 1 fim 2 area. 
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Supplementary Figure 1 : Coalescence of superconducting islands (a) SC island distributions showing individual island size, the number 
of islands of that size, and the fraction of the film area occupied by islands of that size for the substrate of Figs. 1(a), 2(a-e), and 3. SC film 
regions are identified for d(x, y) > d\f T = 0.7 nm on a lxl jim image. Five different <i dep are shown, three of which match those in Figs. 2 
and 3. d dep = 1.07 nm is the 1ST. Vertical dashed lines are plotted to indicate a coherence area, 7r(£/2) 2 with £=15 nm, and the hole-array 
unit cell area, A ce ii =9900 nm 2 . (b) Film area coverage of the largest SC island for each deposition in (a) versus d dep . The 1ST is shown as 
a vertical dashed line at <i dep =1.07 nm. The shaded region indicates a 10% uncertainty in the thickness measurement, (c) The analysis of (a) 
and (b) is repeated for films on all six different NHC substrates and here we plot the deposition thickness at which 50% of the film area is 
covered by a single SC island, d 50 %, compared to the experimental value g£ist for each substrate. The open circle represents the substrate in 
(a) and (b). 



Supplementary Information 

The evolution of the distribution of SC island sizes with thickness provides further support for the islanding picture. Supple- 
mentary Fig. 1(a) presents histograms of the SC island sizes for the experiment of Fig. 1(a) for five thicknesses. For insulating 
films with d dep < 0.89 nm, the island size is relatively uniform and on the order of the estimated coherence area. There as 
many as 71 islands in one bin at 75 nm 2 that together make up 0.7% of the film area at d dep = 0.76 nm. Assuming the islands 
are round (see Fig. 3(a) and (b)), their diameters are ~ 12 nm and 25 nm for d dep = 0.76 and 0.89 nm, respectively. As d dep 
nears g?ist, individual islands coalesce to form larger islands. At the 1ST, d dep = 1.07 nm, a single SC island occupies ~50% 
of the film area. This island grows rapidly when d dep ~ disr as shown in Supp. Fig. 1(b). Denoting d 50 % as the thickness 
at which the dominant island covers 50% of the film area, we find that d 50 % ~ c^ist for films on all six substrates (see Supp. 
Fig. 1(c)). Thus, the coalescence of regions that are thick enough to support SC, as defined by a reference film, is necessary for 
the onset of SC in the NHC films. Finally, because coloration of SC regions in the thickness variation images of Fig. 3 and the 
determination of SC islands in Supp. Fig. 1 are based only on d > d r c ei and do not consider relevant coherence lengths, phase 
coherence between islands is not compulsory and is not universally realized until d = disT- 



